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ABSTRACT 

We examine the ability of a future X-ray observatory, with capabihtics similar to those 
planned for the Constellation-X or X-ray Evolving Universe Spectroscopy (XEUS) mis- 
sions, to constrain dark energy via measurements of the cluster X-ray gas mass frac- 
tion, /gas- We find that /gas measurements for a sample of ^ 500 hot (fcT > 5keV), 
X-ray bright, dynamically relaxed clusters, to a precision of 5 per cent, can be used 
to constrain dark energy with a Dark Energy Task Force (DETF; Albrecht et al. 2006) 
figure of merit of 15 — 40, with the possibility of boosting these values by 40 per cent 
or more by optimizing the redshift distribution of target clusters. Such constraints 
are comparable to those predicted by the DETF for other leading, planned 'Stage 
IV' dark energy experiments. A future /gas experiment will be preceded by a large 
X-ray or SZ survey that will find hot, X-ray luminous clusters out to high redshifts. 
Short 'snapshot' observations with the new X-ray observatory should then be able to 
identify a sample of ~ 500 suitably relaxed systems. The redshift, temperature and 
X-ray luminosity range of interest has already been partially probed by existing X-ray 
cluster surveys which allow reasonable estimates of the fraction of clusters that will 
be suitably relaxed for /gas work to be made; these surveys also show that X-ray flux 
contamination from point sources is likely to be small for the majority of the targets 
of interest. Our analysis uses a Markov Chain Monte Carlo method which fully cap- 
tures the relevant degeneracies between parameters and facilitates the incorporation 
of priors and systematic uncertainties in the analysis. We explore the effects of such 
uncertainties for scenarios ranging from optimistic to pessimistic. We conclude that 
the /gas experiment offers a competitive and complementary approach to the best 
other large, planned dark energy experiments. In particular, the /gas experiment will 
provide tight constraints on the mean matter and dark energy densities, with a peak 
sensitivity for dark energy work at redshifts midway between those of supernovae and 
baryon acoustic oscillation/ weak lensing/cluster number counts experiments. In com- 
bination, these experiments should enable a precise measurement of the evolution of 
dark energy. 

Key words: cosmology: observations - cosmologyxosmological parameters - cosmol- 
ogy: theory ~ x-ray clusters - dark energy 



1 INTRODUCTION 

In the early 1990s, measurements of the baryonic mass frac- 
tion in X-ray luminous galaxy clusters provided compelling 
evidence that we live in a low density Universe. Under the 
assumption that large clusters provide approximately fair 
samples of the matter content of the Universe, X-ray obser- 
vations require that the mean matter density, fim, is signifi- 
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cantly less th an the critical value, with a best-fit value fl 
0.2- 0.3 (e.g.lW hitc fc Frcnk"l99l'; 'FabiaJll99ll; iBriel et aL 



1992 
1995 



White et al. 1993; David ct al. 1995; White & Fabian 
Evrardll 19971: IMohr et al.lll999l: lEttori fc Fabianl|l99' 



19951: lEvrardl 19971: IMohr et al.lll999l:IEttori fc Fabiai 
Roussel et all l20od iGrego et al.1 l200ll: lAllen et al.1 



2004 |200^: lEttori et all l2003l: I Sanderson fc PonmanI 



I999I : 
2002!. 



20031 : 



Lin et al.1 [2OO3I : iLaRogue et al.1 bood ). When combined 
with the expectation from inflation models, later con- 
fir med by Cosmic Microwave Background (CMB) stud- 
ies (|Bennett et al.|[2003l : ISpergel et al]l2003l . and references 
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therein), that the Universe should be close to spatially flat, 
X-ray results on the cluster baryon mass fraction quickly 
lead to the suggestion that the mass-energy density of the 
Universe may be d ominated by a cosmological constant (e.g. 
I White et al.lll993l ). 

The first direct evidence for late-time cosmic accel- 
eration, as would be produced by a sizeable c osmological 
constant, was provided in t he la te 1990s by iRiess et al.l 
ij 19981 ) and lPerlmutter et all (jl999l ) based on measurements 
of the light curves of ty pe la supernovae (SNIa). Since then , 



20031 : |Ricss ct al.1 12004|: 



larger SNIa data s ets (|Knop ct al., ou^., i^^^t,. 

Astier et al] l2006l : iRics s ct al.. l2007l : Iwood-V asev et al l 
20071 : Davis et al.l l2007l ) and an increasingly wide array 
of other, complementary experiments have confirmed 
and improved upon this striking measurement. The com- 
bination of CMB data from the Wilkinso n Mic r owave 
Anisotropy Probe (WMAP) (|Spergel et al.1 |2003| . l2007l : 
iDunklev et al.|[2008h with large scale stru cture (LSS) data 
from t he Sloan Digital Sky Survey (SDSS) (jEisenstein et al] 
I2OO5I : IPercival et al] |2007^ and/or 2d F Galaxy Redshift 
Survey (2dFGRS) (jCole et alj l2005h provides powerful 
evidence for dark energy. The cross-correlation of CMB 
and LSS fiuctuations reveals the effects of dark energy 
on the Integrated Sachs- Wolfe effect (|Scranton et alJ 



I2OO3I : iFosalba et aD l2003l : iRassat et all l2007l l. Mea- 
surements of the amplitude and evoluti on of matter 
fluctu a tions using X-ray ga la xy cl us ters (|Borgani et alJ 
2OOII: iReiprich k. Bohringeil |2002|: lAUen et al.1 12003 " 



Schuecker et al, 
20041: iMantTet al. 



20031 : 



( Gladders et al.1 12007: 

est data jViel et al' 20041: Seliak et al 



Voevodkin fc VikhlininI 1200. 
20081'). optically-s elected clusters 
Rozo et al. I |2007|). Lv man-a for- 
"""" 2OO5I), and weak 



lensing (IVan \yaerbeke et al.1 l2005l: Jjarvis et al.1 l2006l : 



iHoekstra et all I2OO6I : iBeniamin et al. I l2007l ). also provide 
important, powerful confirmation of the new, standard 
cosmological paradigm: a universe in which the main mass 
and energy components are dark matter and dark energy, 
and where dark energy drives the current acceleration. The 
standard model for dark energy remains the cosmological 
constant, which is mathematically equivalent to vacuum 
energy. In principle, however, cosmic acceleration could be 
driven by either dark energy or a m odification to the laws 
of gravity on cosmological scales (see ICopeland et al]|2006l . 
for an extensive review). 

Building o n the early X-ray work , Allen et al.l (|2004f ): 
iRapetti et"S] (120051 ): and lAUen et al] (|2008l ) showed that 
measurements of the evolution of the X-ray gas mass frac- 
tion, /gas, in the largest, dynamically relaxed galaxy clus- 
ters provides a further powerful, complementary approach 
for studying dark energy. As with SNIa data, /gas(z) mea- 
surements probe the redshift-distance relation; whereas the 
peak SNIa luminosity varies as the square of the distance, 
/gas measurements vary as distance, d^ '\ (The distance de- 
pendance derives from the way in which /gas values are 
determined from the observed X-ray temperature and sur- 
face brightness data; lAllen et al.ll200^ In combination with 
the tight constraint on fim provided by the normalization 
of the /gas(z) curve, under the assumption of fair matter 
samples, the /gas(z) data contain sufficient information to 
break the degeneracy between fim and the dark energy equa- 
tion of state, w, in the distance equations. The additional 
combination of /gas and CMB data breaks other important 



degeneracies between parameters in cosmological analyses 
jRapetti et alj|2005l: lAllen et aLlliooj ). 

lAUen et al.l (|2008l ) show that the current constraints on 
dark energy from the /gas experiment are of comparable pre- 
cision to other leading techniques, and are robust under the 
inclusion of conservative systematic allowances, e.g. relaxing 
the requirement for exact hydrostatic equilibrium and al- 
lowing for moderate redshift evolution in the cluster baryon 
fraction. These authors also show that intrinsic, systematic 
scatter remains undetected in the current /gas data, despite 
a weighted mean statistical scatter in the individual distance 
measurements of only ~ 5 per cent ; in contrast, SNIa studies 
iRiess et al.ll2007l : Ijha et al.|[2007l : IWood-Vasev et al.|[2007h 
have established the presence of systematic scatter at the 
~ 7 per cent in distance measurements from the best cur- 
rent SNIa data. 

The key to determining the nature of dark energy is to 
obtain precise measurements of its evolution with redshift, 
z, or sc ale factor, a = 1/(14- 2:). The Dark Energy Task Force 
report (|Albrecht et alj|200d hereafter DETF) presented es- 
timates of the constraints on dark energy parameters that 
should be achievable with a number of future proposed or 
planned dark energy experiments. In particular, the report 
forecasted the ability of these experiments, in combination 
with CMB data from the Planck satellite, to constrain a 
dark energy model of the form w{a) = wo + Wa(l — ti), and 
defined a figure of merit (hereafter FoM) to allow for easy 
comparison of the constraints. In this paper, we use the same 
dark energy parameterization and FoM to quantify the con- 
straining power of future /gas experiments, to be carried out 
with e.g. the Constellation-X or X-ray Evolving Universe 
Spectroscopy (XEUS) missions, in combination with CMB 
data. We show that the /gas experiment is likely to provide 
comparable constraining power to the best other, contempo- 
rary space and ground-based experiments described by the 
DETF. When combined, future CMB, SNIa, baryon acous- 
tic oscillation (BAO), weak lensing, cluster number count 
and /gas experiments should provide precise, accurate con- 
straints on w(z) and allow significant progress in under- 
standing the origin of cosmic acceleration. 

The structure of this paper is as follows: in Section [2] 
we define the dark energy model and the FoM. In Section [3] 
we describe the simulated /gas and CMB data sets. For the 
/gas data, we assume instrument characteristics appropriate 
for the baseline Constellation-X mission. The CMB data set 
approximates that expected from two years of Planck data. 
We also simulate a data set representative of that produced 
by follow-up observations of the Sunyaev-Zel'dovich effect in 
the clusters targeted for the /gas work. Section 3] describes 
the Markov Chain Monte Carlo (MCMC) pipeline and de- 
tails of the analysis method. Our main results are presented 
in Section O Section [S] summarizes our conclusions. 



2 THE DARK ENERGY MODEL AND FOM 

We characterize the evolution of dark energy by its en- 
ergy density in units of the critical density, Side, and its 
equation of state, w. Following the DETF, we parameter- 
ize the evolution of the dark energy equation of state as 
w(a) = wo+ Ma(l - a) (|ChevaUier fc PolarsklbOOll : iLinded 
l2003h for which a cosmological constant has w{a) — —1. In 
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this model, the dimensionless Hubble parameter as a func- 
tion of scale factor has the form 

E{a) = = V^ma-' + nde/(a) + n^a-^ , (1) 

-no 

where 

f{a) = a-3(i+»o+»a)g-3™.{i-«) ^ (2) 

Ho is the present-day value of the Hubble parameter and Sim 
and fik are the mean matter density and curvature density 
in units of the critical density, respectively. 

Using this parameterization, the DETF define a FoM 
that is used to compare the constraining power of difi'erent 
dark energy experiments. Nominally, the FoM scales with 
the inverse of the area enclosed by the 95 per cent confidence 
contour in the wq — u)a plane. However, the DETF showed 
that since there is little correlation in the Wp — lUa plane, 
the area is also proportional to the product of the standard 
deviations o-{wp) x (j(wa), where Wp = w{ap) is the pivot 
value of w{a), i.e., the value of w{a) at which its uncertainty 
is minimized (|Lindeill2006l l. (Note that the standard error 
a{wp) approximately corresponds to the 68.3 per cent uncer- 
tainty in w that would be obtained for a constant— w dark 
energy model). This leads to the definition 

FoM ^ [aiwp) X a{w^)]-\ (3) 

For the DETF Fisher matrix analysis, the la confidence 
region in the Wp — w^, plane forms an ellipse for which 
the semi-axes are the standard deviations of Wp and uia. 
For the more detailed MCMC analysis used here, however, 
we obtain slightly asymmetric probability distributions for 
these parameters in some cases, although to either side 
of the peak probability the distributions can be modelled 
as approximately Gaussian. Therefore, in calculating the 
FoM, we model the la confidence contour in the Wp — Wa. 
plane with a geometrical shape formed by four quarters 
of four different ellipses for which the semi-axes are the 
standard deviations of the Gaussians to either side of the 
peak, namely aup{wp), crdown(wp), aupiwi,), and 
The area of such contour is equivalent to the area of an 
ellipse with semi-axes a{wp) = [aupiwp) + crdown(wp)]/2 
and a{ws) = [crup(wa) -f crdown(«'a)]/2. Thus, we calculate 
our FolVJj as the inverse of the product of the semi-axes 
[a-{wp) X (T(it)a)]~^ which allows a direct comparison with 
the results reported by the DETF. 



3 SIMULATED X-RAY DATA 

3.1 A strategy for future /gas work 

We assume that a future /gas experiment will be carried 
out by an X-ray observatory with capabilities comparable to 

^ To confirm the validity of our definition of the FoM we have ex- 
plicitly measured the area contained by the filled contours in the 
right panel of Figure [2] Dividing this area by both the geometric 
factor TT, which accounts for the conversion between the area of an 
ellipse and a quarter of its circumscript rectangle, and the factor 
2.3, which accounts for the change in the confidence levels from 
two to one degrees of freedom, we successfully match the mea- 
sured area to the value obtained by the product (T{wp) X a{wa). 



Table 1. Baseline X-ray observatory characteristics. 



Band pass 


0.3-10 keV 


Spectral resolution 


E/AE ' 


~ 2400 {m keV) 


Effective area 


15,000 


cm^ (®1.25 keV) 


PSF 


^ 15 arcsec (half power diameter) 


Field of View 


^5x5 


arcmin^ 



those of Constellation-X, as summarized in Table 1. The ma- 
jor improvements of such a mission with respect to current 
X-ray observatories are in collecting area, which is a factor 
~ 100 larger than that provided by the Chandra X-ray Ob- 
servatory, and spectral resolution □ We assume that the /gas 
experiment will be preceded by, and will build upon, forth- 
coming X-ray and/or SZ cluster survey^ that will scan a 
significant fraction of the sky and find a large number of hot. 
X-ray luminous, high— z clusters. These surveys will provide 
the initial target lists for the /gas experiment as well as al- 
lowing an array of complementary cosmological tests based 
on t he power spectrum an d mass function of galaxy clusters 
(e.g. lAlbrecht et~alll2006l ). 

From initial surveys of tens of thousands of clusters, 
the ~ 4000 most X-ray luminous (or highest integrated SZ 
fiux) clusters will be identified. The new X-ray observatory 
will then be used to take short snapshot exposures (~ Iks) 
of these clusters, to identify the most apparently dynami- 
cally relaxed syst ems that are most suitable for /gas work 
(| Allen et al.ll2008l ). The selection of relaxed clusters is likely 
to be based primarily on X-ray morphology, but will also uti- 
lize the high spectral resolution capabilities to measure bulk 
gas motions0 The most relaxed clusters will be re-observed 
with deeper exposures to measure the gas mass fraction to 
the required level of precision. 

Current studies of the Massive Cluster Survey (MACS) 
jEbeling et al.ll200ll . l2007h show that at redshifts z < 0.5 
appro ximately 1/4 clus ters are sufficiently relaxed for /gas 
work (|Allen et al.ll200^ ). We (conservatively) calculate pre- 
dicted cosmological constraints for two separate /gas data 
sets, containing either ~ 500 or 250 relaxed clusters. That 
is, we assume that only approximately 1/8 or 1/16 of the 
4000 hottest, most X-ray luminous clusters detected in a 
future survey will be suitable for use in the /gas experiment. 

For the 500— cluster sample, we assume an average ex- 
posure time per cluster of ~ 20ks. For the 250— cluster sam- 
ple, the typical exposure is ~ 40ks. In both cases, the to- 

For details on planned X-ray observato- 
ries see |http: / / constellation.gsfc.nasa.gov/ 1 and 

|http: / /www.rss d.esa.int/index.php?project=XEUS| 

Forthcoming X-ray survey missions 

include Spectrum- RG/cROSITA; see 

|http: / /www.mpe-garching.mpg.de/projects. html#erosita| 

and |http: / /www.mpe-garching.mpg.de/erosita/MDD-6.pdf[ 
Several large-area SZ surveys are already un derway, in- 
cludi ng the South Pole Telescope (SPT) (e.g. iRuhl et al.l 
I2OO4I . see |http://spt.uchicago.edu/[ l , and the Atacama 
Cosmology Telescope (ACT) (e. g. ISeheal et all l2007t see 
'http: / / wwwphy.princeton.edu/act / ) . 

■* The snapshot observations will also bo of great benefit for a 
range of ancillary cluster science. 
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tal time required to complete the /gas observations will be 
< 15Ms. For the assumed instrument characteristics, we ex- 
pect statistical uncertainties in the /gas measurements re- 
sulting from 20ks exposures of ~ 5 per cent, which corre- 
sponds to ~ 3.3 per cent in distance. For typical exposures 
of 40ks, we expect to measure /gas to ~ 3.5 per cent or 
distance to ~ 2.3 per cent. In Section [5] we show that the 
constraints on dark energy from both the 500 or 250-cluster 
sample are comparable. We adopt the 500— cluster sample 
with 5 per cent /gas measurement uncertainties as our de- 
fault data set. 



3.2 The simulated /gas data set 

3.2.1 The luminosity function of clusters 

To simulate the /gas data set, we first need to predict the 
redshift distribution of clusters. We assume an X-ray fiux- 
limited cluster survey similar to that expected to be pro- 
duced by the Spectrum- RG/eROSITA mission, with a fiux 
limit of Fiim = 3.3 x 10"" erg cm"^ s^Mn the 0.1 - 2.4 keV 
band and a uniform sky coverage of /sky = 0.5. We calculate 
the number of clus ters expected to be observed, M, in each 
redshift bin, zi, as (jMantz et al.ll2008l l 



iVi(^i) = 



dn{M, z) 
dM * 



MM. 



where 



/ dL' 

Jo JL,: 



dLP{L'\M)p{L\L'). 



(4) 



(5) 



Here, V is the comoving volume, n(M, z) is the comoving 
number density of halos with a mass less than AI at redshift 
z, L' is the intrinsic luminosity of a galaxy cluster associated 
with a halo of mass M, and L is its luminosity inferred 
from observations. P{L'\M) is the probability for a cluster 
of mass M to have an intrinsic luminosity L'; p{L\L') is the 
probability for a cluster with intrinsic luminosity L' to be 
observed with luminosity L; and LumC^) is the luminosity 
limit function. We calc ulate the co moving volume element 
per redshift interval as (|Hogg|ll999l ) 



dV 



■4-Kf, 



sky 



C {l + zfdAjzf 

Ho 



E{z) 



(6) 



where c is the speed of light, and cJa the angular d i ameter 
distance. Using N-body simulations, Ijenkins et al.l (|200ll ) 
obtained the following fitting formula for the mass function 
of dark matter halos: 



(7) 



dlncr-i M 

where p is the comoving mean matter density of the Universe 
and A, B and e are fitted parameters. Here a'^{M,z) is the 
variance of the linearly evolved density field, smoothed by 
a spherical top-hat filter, W{k;M). In Fourier-space repre- 
sentation, 



a{M,z) 



27r2 



k'^P{k)W^{k\M)Ak, 



(8) 



where k is the wave number, P{k) is the power spectrum 
of the linear density field extrapolated to redshift zero and 
D{z) is the growth factor of linear perturbations normalized 
to be 1 when z = 0. We calculate the power spectrum using 



the GAME code (|Lewis et al.ll2000lf l. For halo findir i g algo - 
rithms tied to the mean mass density, I Jenkins et al.l ( 200 ll ) 
showed that the values oi A, B and e are almost invariant un- 
der both a broad range of cosmologies and redshift. However, 
these authors also showed that these parameters depend 
on the cluster findi ng algorithm. Here, we use A = 0.316, 
B = 0.67, e = 3.82 (|jenkins et al.ll2001^ . which are appro- 
riate for the spherical overdensity al gorithm SO(k = 324) 
Davis et al ] | 19851 : iLacev fc Colli 1994h , where k is the mean 
overdensity of the halo with respect to the mean matter den- 
sity of the Universe. 

In equation Q we hav e a log-normal probability distri- 
bution (|Mantz et al.ll2008l 'l 



P{L'\M) 



JlogioL'-logioi'(M)]=^/2o 



L' ln(10)V2^cr 



(9) 



where L'{M) is the best fit luminosity for a given mass M, 
and a is it s scatter, determined from t he mass-luminosity 
data set of iReiprich fc Bohringeij (1200 2| ) using the relation 



logi 



M E{z) 

h72 Mq 



ix(0.1 -2.4keV) 



1044/1-32 gj.gg- 



^E{z 



,(10) 



for which a = 0.67 and A = l ogm \Mo/ihZ . }M(^) ] = 14.49, 
and a = 0.12, as obtained bv lMantz et al.l (|2008h . In equa- 
tion ((5| we also have a Gaussian probability distribution 



p{L\L') 



JL-L'f/2af 



27ro"i 



(11) 



with standard deviation ai = {anpi^/nph)L. Here riph is the 
number of photons detected from a cluster in the survey and 
""nph = ywph is the associated Poisson error. We assume that 
at the fiux limit of the survey, -Fum, or equivalently at the 
luminosity limit Li — L(_Fiim,z), the number of photons is 

/"■ph, lim 



^ph.lim ' 



20, 



im, Z) 

Using this, we have ai 



(\/W^ 



3.2.2 Temperature selection 



In or der to minim i ze sy stematic scatter in the /gas experi- 
ment, lAUen et al.l (|2008l ) restrict their analysis to dynami- 
cally relaxed clusters with mean gas mass-weighted temper- 
atures measured within 7-2500 fcT25oo > 5 keV. We impose 
the same temperature cut in this analysis, calculating the lu- 
minosity limit Li t hat corresponds to this temperature limit 
from the relation (|Brvan fc Normanlll998l ) 



log 



Lx(0.1 - 2.4 keV) 



^ A + Blog 



kn \ 
kovj 



(12) 



1044/if|ergs-i£;(2 

where Te is the emission weighted X-ray temperature 0. Fit- 
ting the above re lation (|12[l to the X-ray lumin osity and tem- 
perature data of IReiprich fc Bohringeij (|2002l 'l using the lin- 



^ http://camb.info/ 

^ r25oo is the radius within which the mean density is 2500 times 
the critical density of the Universe at the redshift of the cluster. 
^ Te scales with T2500 as fcTe ~ kT25oo/v with r} ~ 1.1 — 1.2, based 
on MACS clusters spanning the redshift range 0.3 < z < 0.7. 
Beyond redshift 0.7 the value of r) slowly decreases towards ^ 1. 
To be conservative, however, we ignore the difference between Te 
and T2500, i.e. we assume r; = 1, at all redshifts. 
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Figure 1. The redshift distribution (solid curve) of clusters above 
the Spectrum- RG/cROSITA X-ray flux limit with temperatures 
^T2500 > 5keV. A sky coverage of 50 per cent is assumed. This 
redshift distribution has been used to generate the mock /gas 
data set. Assuming that ~ 1/8 of such clusters will be suffi- 
ciently relaxed for /gas work, we obtain a final sample of ~ 500 
clusters. For comparison purposes, we also show (dashed curve) 
the redshift distribution for the case of a fixed luminosity limit 
Lx(0.1 - 2.4keV) > 3.35 x 10''*/if(f ergs"! (no temperature cut) 
which gives a similar total number of clusters. The latter distri- 
bution has more high— z clusters. 



ear re gression BCES(Y|X) algorithm of lAkritas fc Bershadvl 
(|l996l ). we obtain A = -1.46 ± 0.09 and B = 2.50 ±0.1313 
The limiting luminosity in equation ((5)l is then 

^2l 



Liimiz) = min[Li, AnFnn-idi^iz) 



(13) 



with the appropriate if- correction applied in calculating 
the Fiim values. 



3.2.3 The redshift distribution o/ /gas clusters 

Table [2] summarizes the parameters describing our fiducial 
cosmology. For this cosmology, we have calculated the red- 
shift distribution of galaxy clusters over the range < 2 < 2. 
Our fiducial cosmology approximately matches that used 
by the DETF, but includes updated values for Tig and r 
to bet ter match the WMAP three-year and fi ve-year re- 
sults (|Spergel et al.1 l2007l : iDunklev etHI 120081 ). We also 
adopt a lower value for erg = 0.8, consistent with both the 
WMAP three-year and five-year results and the results of 
iMantz et al.l (|2008l ) from measurements of the X-ray lumi- 
nosity function of galaxy clusters within z < 0.7. 

Figure [T] shows the redshift distribution (solid line) for 
clusters detected above the Spectrum-RG/eROSITA X-ray 



flux limit with mass-weighted temperatures kT25oo > 5keV. 
A sky coverage, /sky = 0.5 is assumed. Approximately 5000 
clusters meet these criteria from which, following our ob- 
serving strategy, 4000 will be observed by short snapshots. 
Assuming that ~ 1/8 of these clusters will al so meet the 
relaxation criteria based on X-ray morphology (|Allen et al.l 
I2OO8I : iMillion fc Allenlliooi ). a sample of ~ 500 hot. X-ray 
luminous, dynamically relaxed clusters can be defined. Tak- 
ing snapshot observations of the available ~ 5000 clusters 
instead of 4000, and assuming that ~ 1/8 of these clusters 
are relaxed, we will obtain a sample of ~ 625 /gas targets. 
This allows us to either use a larger sample of clusters, as- 
sume an even more conservative ratio of relaxed clusters, or 
select a different redshift distribution for the /gas sample of 
~ 500 clusters. In Section [5.31 we discuss the latter case. 

For comparison purposes. Figure [T] also shows (dashed 
curve) the redshift distribution for the case of a luminosity 
limit of Li > 3.35 x lO^'^/i^ erg in the 0.1 - 2.4 keV band 
(dashed line; no temperature cut is imposed). The effect of 
the X-ray flux limit on the distribution is evident towards 
the highest redshifts {z ~ 1.5) in this case. 

It is clear from that figure that the temperature and lu- 
minosity cuts lead to different redshift distributions 0. In the 
case of the temperature cut (solid line), the redshift distri- 
bution peaks around z ~ 0.65 and relatively few clusters are 
found at z > 1.5. For the case of the luminosity cut (dashed 
line), the distribution peaks around 2; ~ 1, and has many 
more clusters in the redshift range 1 < 2 < 2. It is impor- 
tant to note, however, that a redshift distribution weighted 
towards higher redshifts does not necessarily imply tighter 
constraints on dark energy. For the DETF FoM criterion, 
constraints around the pivot redshift are important; for the 
/gas experiment Zp ~ 0.25 (see Figure (5)1. In Section [5.31 we 
further discuss the effect that using different redshift distri- 
butions has on the dark energy constraints. 

We generate mock /gas measurements for 500 clusters 
with the redshift distribution appropriate for the case of the 
temperature cut [solid curve. Figure [T] in acco rdance with 
the s election criteria used for current /gas work (|Allen et al.l 
I2OO8I )]. For each cluster, we assign a statistical error in the 
/gas measurements of ~ 5 per cent. We have also generated 
a set of mock measurements for the case of 250 clusters 
observed with /gas measurements accurate to 3.5 per cent. 
This latter data set is used to study the impact on the dark 
energy constraints in the case that the fraction of suitably 
relaxed clusters is less than 1/8 at high redshifts. 

We stress that the predicted redshift distribution, which 
peaks around 2 ~ 0.65 in the case of the temperature cut, 
has already been probed, at least partially, over the lumi- 
nosit y and temperature r ange of interest, by the MACS sur- 
vey (jEbeling et al.l 120011 ): MACS covers the redshift range 



0.3 < 2 < 0.7 to a fiux limit of Fum = 10"^^ ergcm"^ s"^ in 
the 0.1 - 2.4keV band. For MACS, approximately 1/4 
clusters are found to be sufficiently relaxed for /gas work 
(|Allen et al.ll2008l ). Therefore, our assumption that ~ 1/8 
clusters detected in a future X-ray survey and meeting the 
X-ray flux and temperature criteria will be suitably relaxed. 



* We exclude objects from the lReiprich &: Bohringeil ||2002| ) sam- 
ple for which the temper ature was est i mated from the luminosity- 
temperature relation of iMarkevitchI l ll998l) rather than directly 
measured. This leaves 88 data points in total. 



^ The redshift distributions shown in Figure[T]are sensitive to the 
mass-observable relation obtained by iMantz et al.l (|2008ft using 
current data. See that work for details. 
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Table 2. Parameter values of our fiducial cosmology, which is a 
fiat ACDM cosmology. 



wo = -1 


= 0.27 


Ho 


= 72kms-i Mpc-i 


Wa =0 


Ci, = 0.046 


ris 


= 0.95 


Qae = 0.73 


0-8 = 0.8 


A, 


= 2.3 X 10-9 


Hfc = 


b = 0.82 


T = 


= 0.09 



appears reasonable. Moreover, as discussed in Section 5.1, 
for the case of the 250-cluster sample (i.e. assuming that 
only ~ 1/16 clusters are relaxed) and using a similar total 
observing time to obtain individual /gas measurements to 
~ 3.5 per cent accuracy, we obtain very similar dark energy 
constraints (see Table |4]). 

A final important point regards contaminating point 
sources: for MACS clusters, the fraction of the measured 
0.1 — 2.4 keV X-ray flux arising from contaminating point 
sources is small, typically of order a per cent (Mantz et al. 
2008; this is also the case for the hottest, kTc > 5keV, re- 
laxed clusters at lower redshifts.) Therefore, we do not ex- 
pect our target clusters, which have comparable X-ray tem- 
peratures and luminosities, to be severely affected by con- 
taminating point sources, especially at z < 1. This alleviates 
the instrumental requirements on the point spread function. 
An instrument with capabilities similar to the baseline char- 
acteristics listed in Table 1 should be capable of making 
significant strides in dark energy work. 

3.3 Follow-up SZ observations 

The thermal SZ effect is a modification to the CMB spec- 
trum caused by Compton scattering of CMB photons by 
hot electrons in the intracluster medium. The SZ fiux mea- 
sured at radio or sub-mm wavelengths can be expressed in 
terms of the Compton parameter. For a given cosmol- 
ogy, the y— parameter can also be predicted from the same 
X-ray data used to determine the /gas measurements, be- 
ing proportional to the integral along the line-of-sight of the 
product of electron density and temperature, J ricTcdl. 

We have examined the additional cosmological con- 
straining power that can be achieved with follow-up 
radio/sub-mm SZ observations of our sample of 500 clus- 
ters, assuming direct SZ fiux measurements accurate to 2 
or 5 per cent (a level of accuracy that should be straight- 
forward for SZ detector technology ava i lable at the time 
of the experiment; see iMuchovei et al. I l2007l . and refer- 
ences therein). The statistical uncertainties in the predicted 
Compton J/— parameters will be comparable to those asso- 
ciated with the /gas measurements: ~ 5 per cent for the 
500-cluster sample. We generate our predicted j/— parameter 
data set for the redshift distribution shown in Fig. 1 (solid 
curve). 

3.4 Mock CMB data sets 

We have used the CAME code (|Lewis et al.ll200Ct l to gen- 
erate auto and cross temperature and polarization angular 
power spectra, Cj^^, C^^ and Cf^, for the fiduci a,l, flat 
ACDM cosmology described in Table [2] We follow iLewid 



(|2005l ) and iLewis et al.l (|2006h and assume that the tem- 
perature, T, and polarization i?— fields are Gaussian and 
isotropic. We also assume that the polarization B— field is 
negligible. 

Having Ci''^ , C™ a nd C^^, we add a simple, isotropic 
noise power spe ctrum (|Coorav et al.l l2000l : I Lewis! l2005l : 
iLewis et allbOOel ) 

iV,(^',l.)=M(^)^.^.e''M''('+^)/«'■^^ (14) 

where J\fi{iy) = [b{p)a{iy)f, b{u) = VShPlo [v) is the beam 
full width at half maximum (FWHM) measured in radi- 
ans, and (j{u) = {AT/Tf is the root mean square noise 
per beam-sized pixel. Assuming uncorrelated n oise in the E 
and T fields, the covariance over realizations is (|Lewis et al.l 
l200d) 

C l - 1^ ^TE ^EE _^ ^EE j ■ (15) 

For the channel ly = 143 GHz , M^'^ M = = 
2 X 10"* ^iK^ and b{u) = 7.1 arcmin (|Lewis et al.1 120061 ). 
These values correspond to a{v)'^'^ — 6.97 p K and cr(iy)^^ = 
9.68 /iK, which is roughly the sensitivity expected for 
PLANCyj af ter ~ 2 years (14 months) of a full sky survey 
(iPlancklHool . 

We consider two different scenarios relating to fore- 
ground CMB polarization contamination. Firstly, we ex- 
amine the idealized case where such contamination can be 
neglected (iBond et al.ll2004l : iLewidbOOsI : iLewis et al]|2006l : 
|Planckll2006h . Secondly, we consider a more conservative sce- 
nario where ~ 20 per cent of the sky is irretrievably contam- 
inated by Galactic emission, leaving ~ 80 per cent that can 
be modelled as approximately foreg round-free. For the sec- 
ond scenario, iTegmark et al.l (|2000l ) forecast that Planck 
will be able to determine the optical depth to reionization 
to a precision of o"(r) ~ 0.01, as compar ed to cr(r) ^ 0.005 
for the idealized , foreground- free case IjBond et al. 
iLewis et al.]|2006l : |Planckll2006t ). To account for the effects of 
polarization contamination, the DETF discarded polariza- 
tion information for multipoles / < 30 and imposed a prior 
on r to obtain o"(r) = 0.01. For our analysis in the case of 
polarization contamination, we also artificially weaken the 
constraints on r to a precision of cr(r) ~ 0.01 by enlarging 
by an order of magnitude the noise at low multipoles Z < 30 
in the polarization data. 

For both scenarios, we use only the data from mul- 
tipoles 2 ^ Z ^ 2000. For simplicity, we adopt the zero- 
contamination scenario as our default CMB data set. 



4 DATA ANALYSIS METHOD 

4.1 Markov Chain Monte Carlo (MCMC) code 

Given the dark energy model described in Section [5] and the 
simulated /gas and CMB data sets described in Section |3l 
we use the Metropolis Markov Chain Monte Carlo (MCMC) 
algor ithm implemented in the COSMOMC EB (| Lewis fc Bridle! 
|2002| ) package to examine posterior parameter distributions. 

http:/ /www.rssd.csa.int/indcx.php?project=Planck 
http:/ /cosmologist.info/cosmomc/ 
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We use a modified version of the CAMB (|Lewis et al.ll200ol ) 
code to calculate CMB power spectra; this accounts for the 
effects of dark energy perturbations for ev olvinR dark en- 
ergy equations of state jRapetti et alJlioOSl l (see Section|4]5] 
for details). Our modified version of the COSMOMC code 
also incorp o rates the /gas analysis method describ ed by 
lAUen et alj (|2008l ) (see also lRapetti et al.ll2005l . |2007| ). 

Our choice to forecast parameter constraints using a 
full MCMC analysis has some advantages over the more 
widely used Fisher matrix formalisrn (see discussions in 
IPerotto et ai] l2006l : iLewis et al] l2006h . Firstly the shape 
of th e mean log likelihood [see equation (|22|l ] ([Lewis et alJ 
I2OO6I ) in the MCMC analysis encapsulates all of the relevant 
degeneracies between parameters, which is crucial for non- 
Gaussian distributions. Secondly, the fact that our forecasts 
are made using the s ame COSMOMC a nalysis code used to 
analyze current data (|Allen et al.ll2008l ) ensures consistency 
between present and future constraints. Finally, the MCMC 
method allows us to easily and efficiently introduce priors 
and allowances and thereby study the effects of systematic 
uncertainties. 



4.2 X-ray gas mass fraction analysis 

4.2.1 The /gas method 

The X-ray gas mass fraction, /gas, is defined as the ratio of 
the X-ray emitting gas mass to the total mass of a clus- 
ter. This quantity can be determined from the observed 
X-ray surface brightness and the deprojected, spectrally- 
determined gas temperature profile, under the assumptions 
of spherical symmetry and hydrostatic equilibrium. To en- 
sure that these assumptions are as accurate as possible, it 
is essential to limit the /gas analysis to the hottest, most 
X-ray luminous, dynamically relaxed clusters available [Sec- 
tion [ST] for a det ailed discuss i on of the method and current 
measurements see lAUen et al.l (j2008h and references ther ein. 1 
In order to study dark energy, lAllen et all (|2008l ) use 
/gas measurements for a sample of 42 hot (fcT2500 > 5keV), 
X-ray luminous, dynamically relaxed clusters. The /gas mea- 
surements are made within an angle ^2500*^ for each cluster, 
corresponding to r25oo for a reference flat ACDM cosmology 
(with Qm = 0.3 and Ho = 70kms^^ Mpc~^). The /gas mea- 

are related to 



surements in the reference cosmology /g, 
the true values /gas° as 



rACDM / 
/gas V^i 



/ /gas 



ACDMx 



jACDM 
"A 
Jtrue 

"a 



3/2 



(16) 



Non-radiative hyd r odynamical s i mulat ions ijEke et al ] | 19981 : 
iNagai et all l2007l : ICrain et al.l bOOTh suggest that /^^.^ 
is likely to be a pproximately constant in redshift. Thus 
(| Allen et al.ll2008l l. 



true ^ ntrue\ 

(Z; (72500 J — 



lO.5 



bo 



1 + so 



(17) 



wher e sq = O.mh^" (|Lin fc Mohd 12004 IConzalez et~al] 
I2OO7I PI is the observed ratio of the mass in stars (both in 
galaxies and intracluster light) to the X-ray emitting gas 



12 iLin fc Moh^ bOO-j l measured the mass in stars in galaxies, and 
included the intracluster light only as a model. iGonzalez et al.l 



mass, and bo = 0.82 (|Eke et al.lll998l ) is the depletion factor 
for the baryon fraction in clusters with respect to the cosmic 

mean value. 

As discussed bv lAUen" et al.1 (|2008l ). an angular correc- 
tion factor is also required to account for the fact that 
/gar (2; e^soo) needs not be exactly equal to /*™(2; O^g^). 
Observations of large, relaxed clusters show that for the ra- 
dial range of interest, 0.7 < r/r25oo < 1-2, the /gas(r) pro- 
files can be fit by a shallow power-law model with slope 



77 = 0.214 ± 0.02213- Thus, we have 



/true/ /lACDM \ rtruc / /jtruc n 
gas \Z',O2500 J — /gas (Z',02500) 

where 6*2500 = '"2500 /^a, and 



/lACDM 
''2500 



fltruc 
'2500 



nACDM 

'2500 
Zltruc 
''2500 



V[H(z)dA(2)l'^CDM 



(18) 



(19) 



This correction factor is small and can be neglected for 
most analyses of current data, although its inc lusion leads 
to sli ghtly tighter constraints on dark energy (|Allen et al.l 
|2008| ). However, for future experiments of the precision be- 
ing considered here, the inclusion of the angular correction 
term becomes important. 



4.2.2 Allowances for systematic uncertainties 

Following lAllen et al.l (|2008l '). we modify equation (|17|l to 
account for systematic uncertainties in the /gas analysis: 



/true / /itruc 
gas 172500 



— K I 



b{z) 



l + s{z) 



(20) 



Here 7 allows for departures from the assumption of hy- 
drostatic equilibrium, due to non-thermal pressure sup- 
port; K is a normalization uncertainty relating to instru- 
mental calibration and certain modelling issues; b{z) = 
60(1 -I- QbZ + /3bZ^) accounts for uncertainties in the cluster 
depletion factor, both in the normalization, &o, and possi- 
ble linear, ah, and quadratic, /3b, evolution with redshiflF^: 
s{z) = so{l + a^z + Pbz'^) accounts for uncertainties in the 

stellar mass fractionl^ 

Using hydrodynamic N-body simulations iNagai et al.l 

i2003) show that for measurements at r25oo in large, re- 
laxed clusters, non-thermal pressure support is unlikely to 
exceed 8 per cent. Furthermore, if, as suggested by some 
current X-ray data JPabian et all l2003l : iFabian et all l2005l : 
iRevnolds et al.l 120051 ). the gas viscosity is higher than that 
included in current simulations, then non-thermal pressure 
support could be even lower. Based on these findings, we 



ll2007h measured the mass in stars in both galaxies and intraclus- 
ter light. For the largest clusters both works find simular results. 

Note that even using two very different reference cosmologies 
such as SCDM (Om = 1, Ho = SO kms'^ Mpc-^) an d ACDM 
(Hm = 0.3, Hq = 70kms-l Mpc~M. lAllen et al.l l l2008l) obtained 
similar values for rj around ^2500 . 

^■^ Note that the allowances on ctb and /3b can also be assumed to 
encompass the combined uncertainties in the redshift evolution of 
7, K and 6, which have the sam e effect on equatio n I I20I I. 

Working with current data, lAUen et al.l bOOSi ) use only the 
linear order of the redshift expansions for their systematic al- 
lowances i.e. csb aud Qs- 
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Table 3. Systematic allowances incorporated in the /gas and XSZ experiments. 



Cluster 



Parameter 



Allowance (optimistic/standard/pessimistic) Type 



Calibration/Modelling K 

Non-thermal pressure 7 

Gas depletion: norm. 60 
Gas depletion; evol. (linear) 

Gas depletion: evol. (quadratic) /3b 

Stellar mass: norm. sg 

Stellar mass: evol. (linear) cng 

Stellar mass: evol. (quadratic) /3s 



1.0 ±0.02/ ±0.05/ ±0.10 Gaussian 

0.96 < 7 < 1.04/0.92 < 7 < 1.08 uniform 

0.82 X (1 ± 0.02/ ± 0.05/ ± 0.10) uniform 

±0.02/ ±0.05/ ±0.10 uniform 

±0.02/ ±0.05/ ±0.10 uniform 

0.16 X (1 ± 0.02/ ± 0.05/ ± 0.10) Gaussian 

±0.02/ ±0.05/ ±0.10 uniform 

±0.02/ ±0.05/ ±0.10 uniform 



XSZ EXPERIMENT 

Calibration/Modelling 
evolution (linear) 



ko 

Ok 



1.0 ±0.02/ ±0.05 
±0.02/ ±0.05/ ±0.10 



Gaussian 
uniform 



adopt by default a uniform prior such that non-thermal pres- 
sure support lies in the range — 8 per cent (although a more 
pessimistic range of — 16 per cent is also considered). Since 
the use of an asymmetric prior would bias the analysis, lev- 
ering Qra above the fiducial value, we employ an equivalent, 
rescaled symmetric prior such that 1— (a/2) < 7 < l + (a/2), 
where a = \1 - 1.08|/1.04. 

The depletion parameter, 60, reflects the thermody- 
namic history of the X-ray emitting cluster gas. Using non- 
radiative simulations of hot, massive clusters of comparable 
size to the real c lusters to be used in the /ga s experiment, 
Eke etal.1 (Il998l) (see also lAUen et al.1 |2004| : iNagai et al] 



20071 : ICrain et al.1 12007| ) obtained 60 = 0.82 ± 0.03 at the 
radius of the measurements r25oo (~ 0.25rvir) and found no 
evidence for redshift evolution: Oh = 0.00 ± 0.03 for mea- 
surements made at r ~ O.Sr^j r, spanning the red shift range 
< z < 1. As discussed by lAUen et al.l (|2008l ). however, 
systematic uncertainties are associated with current predic- 
tions for b{z), due to limitations in the accuracy of the phys- 
ical approximations employed in the simulations. Estimat- 
ing the residual uncertainties in the prediction of b{z) that 
will be appropriate at the time of a future /gas data set 
(~ 2015 — 2020) is difficult. We have chosen to use a range 
of values that extend from optimistic to pessimistic scenarios 
(see Table EJ. 

Current optical and near infrared data for low- 



to- intermediate redshift clusters give sp 



jFukugita et al.1 Il998l : iLin fc Mohif 12004 IConzalez et"al] 
I2OO7I ). Although, at present, the constraints on s{z) for clus- 
ters at 2 > 0.5 are sparse, we expect the form of s{z) to be 
relatively well understood by the time of the /gas experi- 
ment. 

In order to keep the interpretation of our results simple, 
we present results for three sets of systematic allowances: 
for the parameters, if, feo, ctb, /3b, so, Os, /3s, we employ al- 
lowances of either ±2 per cent (optimistic), ±5 per cent 
(standard), or ±10 per cent (pessimistic). In all cases, we 
employ uniform priors with the exception of K and sq, for 
which Gaussian priors are more appropriate and therefore 
used. As noted above, a uniform allowance of ±4 per cent on 
7, is included by default, although the effects of doubling the 



uncertainty in this parameter are also examined. We stress 
that whether 7 = 1 precisely, or Ob, Os etc are precisely zero, 
is not of primary importance to a future analysis: if known, 
the exact values can be incorporated into the default model. 
It is the uncertainties in the values that affect the accuracy 
and precision of the dark energy constraints. 

4.3 Analysis of the SZ data: the XSZ experiment 

For the true, underlying cosmology, the measurement of 
the Compton ?y— param eter from both the X-ray and SZ 
data should match (e.g . Molnar et al.|[200^ : ISchmidt et al.l 
I2OO4I : iBonamente et al]l2006l ) . For a given cosmology the y- 
parameter predicted by X-ray data depends on the square 
root of the angular diameter distance to the cluster, d^^, 
whereas the observed SZ fiux at radio or sub-mm wave- 
lengths is independent of the cosmology assumed. Combin- 
ing the y— parameter results, we can measure the distances 
to the clusters as a function of redshift and, therefore, con- 
strain dark energy. 



ACDM SZobs , / X 

y =y k{z) 



Jtruc 

"a 



1/2 



(21) 



Here y^'-^^''^ is the X-ray measurement of the y-parameter 
for the reference cosmology and y^^°^^ is the radio/sub- 
mm observation^ Following a similar approach to that 
adopted with the /gas data, we incorporate systematic al- 
lowances into equation (|2ip : k{z) = ko{l + akz) accounts 
for the combined systematic uncertainties in the X-ray and 
SZ data y— parameter measurements due to calibration, ge- 
ometric effects, gas clumping, etc., and their evolution. We 
employ Gaussian priors on ko of size 2 (optimistic) or 5 
(standard/pessimistic) per cent and uniform priors on ak of 
size 2 (optimistic), 5 (standard) or 10 (pessimistic) per cent. 

We note that the best clusters to observe for the XSZ 
experiment are the same systems used for the /gas experi- 
ment: the largest, most dynamically relaxed clusters. These 

The limitations of existing SZ data have to date restricted the 
XSZ experiment to me asurements of the Hubble constant (e.g. 
IBonamente et al1l2006l . and references therein). 
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are the clusters for which the SZ signals are strongest and 
for which systematic uncertainties associated with geometry 
and thermodynamic structure are minimized. Note also that 
no additional X-ray observations are required to carry out 
the XSZ experiment, once the /gas data are in hand. 

4.4 Incorporating the CMB data 

In addition to the dark energy model parameters and the 
/gas parameters discussed in Section [4. 2 1 we vary the follow- 
ing eight CMB-related parameters in the MCMC analysis: 
the mean physical baryon density, Qhh^', the mean physical 
cold dark matter density, Odm^'^; the (approximate) ratio 
of the sound h orizon at last scatterin g to the angular diam- 
eter distance (|Kosowskv et al.ll2002l ). 9b; the optical depth 
to reionization (assumed to occur in a sharp transition), r; 
the mean curvature density of the Universe, Slk; the scalar 
adiabatic spectral index, ris; and the scalar adiabatic am- 
plitude, ^s, at fc = 0.05 Mpc~^. We employ a uniform prior 
on ln(j4s). The combination of 9^ and In(As) as parameters, 
rather than Ho and ^s, leads to a more Gaussian prob- 
ability density distribution which, in tu rn, aids sampling 
(|Kosowskv et al.ll2002l : iLewis et aLlbOOd ). 

The degeneracies between dark energy model param- 
eters and Ok are of particular importan ce in the analysis 
jRapetti et al.ll2005l : [Clarkson et al.ll2007l )FI For their fore- 
casts, the DETF include Planck priors in their Fisher ma- 
trix analysis, approximating the role of future CMB con- 
straints as well as the degeneracies between the dark energy 
parameters and Slk- Here, we account fully for the degenera- 
cies between parameters, and the complementarity of the 
data sets. 

Given the vector, e, of CMB-related parameters, we 
sample th e exponential of the following mean logarithmic 
likelihood (|Lewisll2005l : iLewis et al.lbOOel ') 

{lnP(e|eo)) = -i[rr(C,(6o)C,(6)-^)+ln|C7,(6)|], (22) 

where eo is the vector formed by the corresponding fidu- 
cial values of Table [2] Note that where the posterior is non- 
Gaussian, the marginalized constraints on individual param- 
eters need not peak exactly at the fiducial values, although 
for all cases considered here the differences are very small. 

4.5 Dark energy clustering 

Fo r a dark energy model with a constant equation of state, 
w. IWeUer fc Lewis! l|2003l ') and iBean fc Dorel (|2004 ) showed 
that dark energy clustering can have a non-negligible im- 
pact on the constraints, driven primarily by the effect of 
such perturbations on the Integrated Sachs- Wolfe (ISW) ef- 
fect. ISpergel et all (|2007l ) showed that accounting for dark 
energy clustering has a large effect on current dark en- 
ergy constraints derived from CMB data alone. Since, for 
constant— w models, combining the CMB data with e.g. 



l'^ lAllen et all ||2004|2008|) and lRapetti et al.l | |2005|) showed tliat 
the combination of /gas plus CMB data allows one to drop both 
the assumption of flatness and the priors on and h that 

would otherwise be required for the /gas analysis. The /gas-|-CMB 
data combination also alleviates other important parameter de- 
generacies, e.g. between Cl^Jh^, ris and t. 



distance measurements from type la supernovae or X-ray 
galaxy clusters leads to tight constraints on w and a re- 
sult consistent with a cosmological constant {w = — 1, 
for which no dark energy clustering occurs), the impor- 
tance of accounting for dark energy perturbat i ons is reduced 
JWeller fc Lewis! !2003! : !Rapetti et al] !2005! : !Spergel et all 
I2OO7I ). However, w hen one considers mo re general models 
in which w evolves, iRapetti et al] (|2005l ) showed that even 
with the best current data combinations, accounting for the 
effects of dark energy perturbations is important: the con- 
straints on w{a) increase by a factor of ~ 2 with respect to 
the case where dark energy clustering is (wrongly) ignored. 

For our analysis, we assume that dark energy is an im- 
perfect fluid where dissipative processes generate entropy 
perturbations. As suggested by quintessence scenarios, we 
assume a constant, general (non-adiabatic) sound speed 
Cg = 1 in the comoving frame of the fluid (denoted by the 
circumflex ' ). This is the only fram e for which the gen - 



eral soun d speed is gaug e invariant (iBean fc Pore! 



2004) 



20041 ) ■ 



Following IWeller fc Le^ (|2003! ) and iBean fc Dori ., . 

!Rapetti et al.! (12005! ) extended the dark energy perturba- 
tion equations to account f or an evolving d ark en ergy equa- 
tion of state, «;(a). As in !Rapetti et al.l (|2005! ) we calcu- 
late the density, S, and v elocity, v, perturbation equations 
(|Ma fc Bertschingei![l995l ) in the synchronous gauge 



-3H(Cs - 



w)S~{l + w){kv + 3B)+£{w) (23) 
kc^S 



1 



(24) 



where both derivatives, denoted by dots, and the Hubble 
parameter, TL, are with respect to conformal time. B = 5a/a 
is the metric perturbation and 5 is the density perturba- 
tion in the comoving frame of the dark energy fluid. The 
density perturbation, 5, can be recast into the CDM co- 
moving frame density and velo city perturbations , S an d 
V, using the relation given by !Kodama fc Sasa^kH l)l984| ). 
5 — 5 + 3H{1 + w)v/k. Using this relation i n equation (|23||. 
we re co ver the pertu r bation equations of !Weller fc Lewis! 
(j2003l ): !Bean fc Dor3 (j2004! ) except that here w depends 
on the scale factor a and introduces a new source term, 
£{w) = SHthv/k, in the density perturbation equation. Note 
that this term depends on the derivative of the equation of 

state, w. Equation (|24p does no t have a new term. 

As discussed by !Vikman! ()2005h : !Caldwell fc Doran! 

(|2005l ) a single classical scalar field cannot evolve from a 
quintessence-like , w > —1, to pha ntom-lik e, w < —1, be- 
havior. However. lOnemh fc Wooda rd (2002;. !2004! 'l proposed 
a single scalar field model where a super-accelerated phase 
{w < —1) of the c osmic expansion can be a chieved via quan- 
tum effects. Later. !Kahva fc Onemli! (|2007!) show e d that this 
model is stable. Alternatively, ^eng et ahl (|2005! ) ; !Guo et aP 
(|2005h :[H^ (j2005! '): IZhao et all (|2005h suggested the so-called 
quintom model: in this model, the effective equation of state 
of two combined scalar fields, one with w > — 1 and the other 
with w < —1, can cross the cosmological-constant bound- 
ary, = — 1, as it evolves in time. Other models that allow 
w{a) to cross this boundary have also been proposed. How- 
ever, for practical purposes, using an effective, evolving dark 
energy equation of state produces a well-known divergence 
in equati on (I24II wh e n w(a) — —1. This dive rgence can be 
avoided (|Huev!!2004l : !CaldweU fc Doran! !2005^ by imposing 



© 2007 HAS, MNRAS OOO.ITlfTel 



10 D. RapetU, S. W. Allen & A. Mantz 




Figure 2. (Left panel) The 95 per cent confidence contours in the uiq — Wg, for the default dark energy model using the optimistic (2 per 
cent; blue, solid contour), standard (5 per cent; dashed contour) and pessimistic (10 per cent; red contour) allowances. (Right panel) The 
68 per cent confidence contours in the w)a — Wp plane for the default dark energy model using the optimistic (2 per cent; purple, solid 
contour), standard (5 per cent; dashed contour) and pessimistic (10 per cent; red contour) allowances. The marginalized Irr confidence 
intervals on Wn and Wp are used to calculate the FoM. The figure confirms that Wn and Wp are not strongly correlated, as assumed in 
the definition of the FoM (Section 2). 



5 = and V — within the logarithmic singularity region, 
10 — —1 + |e|, where e is infinitesimally small. Inaccuracies 
in this approximation have a negligible impact on the re- 
sulting CM B power spectra (|Rapetti et al.l l2005l : IXia et al] 
I2OO6I . I2OO7I ). 

In what follows, we present results for cases where dark 
energy clustering is either accounted for or ignored in the 
CMB analysis (Section 15. 4p . Dark energy clustering does 
not have a significant impact on the /gas analysis. 



5 RESULTS 

As described in Section[2l to enable a direct comparison with 
the predicted dark energy constraints for other planned ex- 
periments, we parameterize our results in terms of the DETF 
FoM. We present results for a fiducial /gas-l-CMB data set, 
incorporating the statistical uncertainties and systematic al- 
lowances described above, and with zero scatter about the 
fiducial curves. The absence of scatter in the simulated sam- 
ples ensures that the peaks of the posterior probability dis- 
tributions occur at the expected values, in the same way 
that the DETF Fisher matrix analysis does. That is, in or- 
der to compare our results with the DETF, we select the 
same realization as the DETF. Note also that the use of a 
zero-scatter realization does not affect the FoM. We have 
explicitly confirmed this, and that using other realizations 
does not have a large impact on the FoM, by comparing the 
FoM for the fiducial realization to a series of Monte Carlo 
simulations, in which appropiate scatter about the fiducial 
/gas curve was included. 



5.1 Constraints on the FoM 

In the first case, we determine constraints for our 'default' 
analysis: this involves /gas data for 500 clusters measured 
to 5 per cent accuracy and CMB data with negligible fore- 
ground contamination. We ignore the effects of dark energy 
clustering (as do the DETF) and allow linear evolution in 
b{z) and s{z). No follow-up SZ data are included. 

The constraints from the default analysis are shown in 
Figures[2j|4l and [5] The left panel of Figure[2]shows the well- 
known degeneracy between wq and Wa. The right panel of 
that figure shows the constraints in the Wa — Wp plane. The 
results from the MCMC analysis confirm that Wp and lUa 
are approximately uncorrelated, which facilitates the simple 
calculation of the FoM—[a{'Wp) x a{wa)]~^ a described in 
Section [2] 

Table U summarizes the results on Q,dc, wq, lUa, 
Wp and the FoM for the default analysis and 2, 5 or 10 per 
cent systematic allowances. Also included in the table are 
results for five further slightly modified, interesting scenar- 
ios: for the case where we include dark energy clustering 
(Section 14. 5|) : for the case where we use the more conser- 
vative CMB data set (Section 13. 4p : for the case where we 
allow quadratic redshift evolution in the gas depletion fac- 
tor, /3b, and the stellar fraction, /3b (Section 14. 2. 2|l : for the 
case we use the 250-cluster /gas data set with 3.5 per cent 
measurement errors; for the case where we double the sys- 
tematic allowance on 7; and for the case where we include 
extra information from the XSZ experiment (Section 14. 3p . 
For each of these scenarios, we list results for 2 per cent (op- 
timistic), 5 per cent (standard) and 10 per cent (pessimistic) 
allowances. 

Interestingly, we see that constraints on dark energy are 
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Table 4. The la uncertainties on the dark energy parameters and FoM. Systematic allowances of 2 per cent (optimistic), 5 per cent 
(standard) or 10 per cent (pessimistic) have been used. Results are presented for the default model and for six other cases, described in 
the text. We obtain a FoM in the range 34 — 43, for the optimistic allowances, in the range 21 — 33, for the standard, and in the range 
15 — 29, for the pessimistic. 



Allowances 


Run 

Model 


a-(Qm) 


Dark energy parameters 
(T(Cdc) o-{wo) d-{wp) 


O-(tOa) 


FoM 

[a{Wp) X CT(«)a)]"^ 


AFoM/FoM 
(percentage) 


2% 


Default 


0.014 


0.011 


0.130 


0.050 


0.52 


38.5 




2% 


DE clustering 


0.014 


0.010 


0.137 


0.048 


0.56 


37.2 


-3.4% 


2% 


CMB conservative 


0.015 


0.011 


0.135 


0.052 


0.57 


33.7 


-12.5% 


2% 


Quadratic 


0.014 


0.010 


0.146 


0.050 


0.56 


35.7 


-7.3% 


2% 


Half sample 


0.014 


0.011 


0.125 


0.049 


0.51 


40.0 


-1-3.9% 


2% 


Double 'y 


0.019 


0.013 


0.129 


0.055 


0.52 


35.0 


-9.1% 


2% 


Adding XSZ 


0.013 


0.010 


0.122 


0.047 


0.50 


42.6 


-1-10.6% 


5% 


Default 


0.022 


0.016 


0.137 


0.056 


0.58 


30.8 




5% 


DE clustering 


0.020 


0.016 


0.152 


0.055 


0.69 


26.4 


-14.3% 


5% 


CMB conservative 


0.023 


0.017 


0.147 


0.058 


0.66 


26.1 


-15.3% 


5% 


Quadratic 


0.022 


0.016 


0.157 


0.079 


0.61 


20.8 


-32.5% 


5% 


Half sample 


0.022 


0.016 


0.132 


0.058 


0.57 


30.2 


-1.9% 


5% 


Double "7 


0.024 


0.018 


0.137 


0.058 


0.58 


29.7 


-3.6% 


5% 


Adding XSZ 


0.021 


0.015 


0.132 


0.053 


0.57 


33.1 


+7.5% 


10% 


Default 


0.033 


0.024 


0.143 


0.064 


0.62 


25.2 




10% 


DE clustering 


0.029 


0.022 


0.164 


0.059 


0.73 


23.2 


-7.9% 


10% 


CMB conservative 


0.038 


0.027 


0.163 


0.069 


0.73 


19.9 


-21.0% 


10% 


Quadratic 


0.033 


0.025 


0.173 


0.106 


0.64 


14.7 


-41.7% 


10% 


Half sample 


0.033 


0.024 


0.137 


0.064 


0.60 


26.0 


+3.2% 


10% 


Double 7 


0.034 


0.025 


0.144 


0.065 


0.61 


25.2 


0.0% 


10% 


Adding XSZ 


0.028 


0.020 


0.134 


0.059 


0.59 


28.7 


-1-13.9% 



similar for most cases of interest (this is also shown graphi- 
cally in Figure [3)). With the optimistic, 2 per cent systematic 
allowances, a FOM in the range 34 — 43 is obtained. For the 
standard, 5 per cent systematic tolerances, the FOM lies 
in the range 21 — 33. Even with the pessimistic 10 per cent 
systematic allowances, we obtain a FoM in the range 15 — 29. 

The final column of Table |4] summarizes the percentage 
differences in the FoM with respect to the default model 
for each case of interest. We see that, using the standard 
or pessimistic (5 and 10 per cent allowances) the greatest 
impact on the FoM occurs by allowing quadratic evolution 
in the systematic allowances; in this case a ~ 30 — 40 per 
cent reduction in the FoM with respect to the default model 
is observed. Using the standard, 5 per cent allowances, we 
see that accounting for dark energy clustering has only a 
small effect (~ 15 per cent); the inclusion of the XSZ data 
only leads to a modest improvement in the FoM (~ 10 per 
cent). Doubling the uncertainty on 7 does not have a ma- 
jor effect on the results, and for the pessimistic scenario, 
becomes negligible. 



5.2 Comparison with DETF results 

Comparing our results on the FoM with th ose reported by 
the D ETF (page 77 of the DETF report; lAlbrecht et~aLl 
120061 ) ■ we find that the /gas experiment has similar dark 
energy constraining power to other leading, future (DETF 
stage IV) ground or space-based experiments. 

Figure |4] shows the 95 per cent confidence constraints 
in the Qdc — Wp plane for the default model and optimistic 
(2 per cent), standard (5 per cent) and pessimistic (10 per 
cent) allowances. The size of this confiden ce region is in- 
verse ly proportional to the FoM. The DETF (|Albrecht et all 
I2OO6I ) present similar figures, with the same axis scaling, for 
the other, future dark energy experiments. The compara- 
ble constraining power and complementary nature of the 
/gas and other experiments can (at least in part) be seen 
by comparing these figures. In particular, the power of the 
/gas(-fPlanck) experiment in constraining fide is evident. 



The /gas experiment described here would fall under the cat- 
egory of stage IV experiments, as defined by the DETF. 
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Figure 3. Posterior probability distributions for Wa for the cases 
described in Table |4] default scenario (black, solid line), with 
dark energy clustering (blue, long dashed line), conservative CMB 
data (red, dotted line), including quadratic evolution allowances 
(green, long dashed line), using the 250-cluster sample with 3.5 
per cent measurement errors (magenta, dot-dashed line) , doubling 
the systematic allowance on 7 (orange, long-short dashed) and 
adding the XSZ experiment (cyan, long dot-dashed line). The 
(optimistic) 2 per cent systematic allowances are used in every 
case. 

-0.7 1 . . . . 1 



-0.8 




§.6 0.65 0.7 _ 0.75 0.8 0.85 
DE 

Figure 4. The 95 per cent confidence contours in the — Wp 
plane for the default dark energy model and optimistic (2 per 
cent; blue, solid contour), standard (5 per cent; dashed contour) 
and pessimistic (10 per cent; red contour) allowances. The axes 
are scaled to cover the same region as the figures presented by 
the DETF. 



5.3 The relevance of the redshift distribution 

To calculate the dark energy results presented in Table |4j 
we employ a redshift distribution of clusters drawn from a 
simulated X-ray lu minosity function (based on the work of 
iMantz et al ] |2008l . as discussed in Section [3.2.1|) for both a 
given fiducial cosmology and a future, planned X-ray clus- 
ter survey. As discussed in Section [3.2.21 we select clusters 
from this distribution using th e same criterion (k T > 5keV) 
than we use for current data (| Allen et al.ll200a ). To apply 
this selection criterion to the simulated data set we use a 
lumin osity-temperature relation o btained from present-day 
data (jReiprich fc Bohringeil I2OO2I ). Note that in this rela- 
tion the temperatures a re emission weight ed instead of mass 
weighted, as they are in Iau en et al.l (|2008l ) , which makes our 
selection even more conservative. 

Finally, we select only relaxed clusters by scaling the 
distribution with a factor 1/8 (or 1/16 for the half sample 
scenario). This factor is especially conservative at low red- 
shifts {z < 0.5), where the MACS survey has shown that 1/4 
clusters are sufficiently relaxed for /gas work. This suggests 
that selecting a distribution with more clusters at low red- 
shifts is a plausible alternative. Furthermore, after all these 
conservative cuts we obtain a total number of clusters suit- 
able for /gas work larger than 500. This offers us additional 
freedom to build an alternative redshift distribution. In this 
section, we use two alternative distributions to assess the 
impact that choosing a particular distribution has on the 
dark energy constraints. 

Within the limits of the analysis described above, we 
design two test distributions such that their peaks are at a 
lower redshift, z ~ 0.5, than that of the original distribution, 
z ~ 0.65. From each test distribution we form an /gas sam- 
ple of ~ 500 clusters, and ~ 5 per cent /gas measurement 
errors per cluster, as we did for the original distribution. The 
first test distribution is approximately gaussian, with a large 
number of clusters at low redshift, a small tail at high red- 
shifts, and very few clusters beyond redshift 1. The second 
test distribution is similar to the original, but shifted to- 
wards lower redshifts. At high redshifts this distribution has 
less clusters than the original but significantly more clusters 
than the first test distribution. 

Using our default dark energy model, and the 2 per cent 
set of systematic allowances, the first test distribution pro- 
vides an increase in the FoM of ~ 10 per cent with respect to 
the original distribution. As mentioned in Section [3.2.31 the 
pivot redshift, Zp, is important for the DETF FoM criterion. 
Figure[5]shows that for the /gas experiment, Zp ~ 0.25. Inter- 
estingly, the second test distribution provides an increase in 
the FoM of ~ 40 per cent, which suggests that having enough 
high redshift clusters is also important for the DETF FoM. 
Note also that to obtain the /gas measurements for each of 
the test samples, we will require a shorter total exposure 
time than that estimated for the original sample (~ 15Ms). 
Thus, using the same exposure time, an even larger increase 
in the FoM should be achievable with such samples. 

These results indicate that a further analysis is required 
to determine the optimal redshift distribution of clusters 
with which to carry out future /gas experiments. Such anal- 
ysis is beyond the scope of this paper, but we will pursue it 
in a forthcoming publication. 
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5.4 The evolution of dark energy 

Figure [5] shows the evolution of the dark energy equation 
of state as a function of scale factor, 'w{a). The pivot scale 
factor, Op, is the scale factor at which we obtain the tight- 
est constraint on w [that measurement being a{'Wp)]. As 
shown in Figure [S] for the /gas experiment we measure a 
pivot scale factor of ftp ~ 0.8, which corresponds to a pivot 
redshift Zp ~ 0.25. Interestingly, these values lie between 
the pivot scale factors/redshifts reported by the DETF for 
SNIa experiments (op ~ 0.93/zp ~ 0.075) and galaxy clus- 
ter number counts, weak lensing and BAO experiments 
(op ~ 0.65/zp ~ 0.54). Pinning down the evolution of w 
over a wide redshift range will be a crucial for unraveling 
the nature of dark energy. Our results argue that a combi- 
nation of /gas, CMB, BAO, SNIa, weak lensing and galaxy 
cluster number count experiments is likely to prove powerful 
in this regard. 



5.5 The CMB data and early dark energy 

The acoustic s cale at last scatte ring, ia, is tightly constrained 
by CMB data (|Page et al.l2003l ') and is highly sensitive to the 
amount of dark energy at recombination, ia changes drasti- 
cally if the early d ark energy de nsity exceeds the matter plus 
radiation density (|Wrightll2007l ). CMB constraints on /a pro- 
vide a strong constraint on d ark energy parameters at early 
times ijPoran fc Lillevll2002h. defining a wel l -known bound- 
ary in the wp — Wb, plane (|Rapetti et al.ll2005l : lUpadhve et al.l 
l2005l : IWrightl2007l ') R At late times, for our experiment, dark 
energy is constrained primarily by the /gas data, with a small 
contribution from the Integrated Sachs- Wolfe (ISW) effect 
in the CMB. 

A simple exercise provides further insight into how the 
CMB data help in constraining dark energy. For this, we 
re-examine the constraints in the wo — Wa plane obtained 
from the /gas-l-CMB data; the 68 and 95 per cent confidence 
contours for the default model with 5 per cent allowances 
are shown (dashed curves) in Figure (6] The combination of 
/gas-l-CMB data provides tight constraints on Qhh^ , ^dmh^ 
and la. (driven primarily by the CMB data), and on h (driven 
by the combination of both data sets). Using these con- 
straints as priors, we examine the constraints in the wq — uia 
plane that can be obtained from the /gas data alone; the 
results are shown as the red, solid curves in Figure (6] We 
see that the priors encompass some of the CMB constrain- 
ing power, in particular in defining the characteristic upper 
boundary in the wq — Ws, plane. However, they do not con- 
tain the full in formation on, e.g., t he covariance of Q,hh^ , 
VlAmh^ and h (|Rapetti et all l2005l : IWrightl [2007h which is 
also important in constraining dark energy at later times. 

We note that the prior on ia provides a tight constraint 
on the curvature. The blue dotted curves in Figure |B] show 



The presence of the boundar 



the 



wo 



^^—^—^—^^^-^ '-y tfte uiQ — Wa plane 

llRapetti et al.ll2005l : lUpadhve et al.ll2005l: IWrieht|[2007l) makes it 
important to consider, as we do here, simulations that account 
fully for measurement uncertainties but which do not scatter 
about the fiducial curve. Otherwise, scatter towards the CMB 
boundary would increase the FoM, and scatter away would de- 
crease it, compHcating the interpretation of results. 




Figure 5. The Icr confidence constraints on the evolution of the 
dark energy equation of state as a function of scale factor w{a). 
Results are shown for the default model (Table [TJ using the opti- 
mistic (2 per cent; shaded, purple region), standard (5 per cent; 
dashed line) and pessimistic (10 per cent; solid, red line) system- 
atic allowances. The tightest constraints on w{a) occur at the 
pivot scale factor, Op ~ 0.8 (zp ~ 0.25). 



the constraints obtained from the /gas data alone, using only 
the priors on flhh^ and h and assuming flatness. 



5.6 The importance of the XSZ experiment 

The XSZ technique provides a complementary and indepen- 
dent experiment to measure dark energy. Although the in- 
clusion of constraints from the XSZ experiment leads to only 
modest formal improvements in the FoM with respect to 
the results for the /gas -f CMB data (Table [l] as can be ex- 
pected given the relatively weak dependence on dark energy 
in equation 21), it is important to note that the XSZ ex- 
periment relies on different assumptions and has different 
systematic uncertainties. In particular, the XSZ experiment 
is independent of assumptions regarding hydrostatic equi- 
librium, the depletion factor, and the stellar mass fraction. 
Thus, the combination of data from the /gas and XSZ tech- 
niques can help to ensure robustness in the results. In prin- 
ciple, the inclusion of XSZ data can also allow some of the 
priors in the /gas experiment to be relaxed. 



6 CONCLUSIONS 

We have examined the ability of a future X-ray observatory, 
with capabilities similar to those planned for Constellation- 
X, to constrain dark energy via the /gas experiment. We find 
that /gas measurements for a sample of 500 hot (fcr2500 ^ 
5keV), X-ray bright, dynamically relaxed clusters, with a 
precision of ~ 5 per cent, can be used to constrain dark en- 
ergy with a FoM of 15—40. These c onstraints are compar able 
to those predicted by the DETF (|Albrecht et all 120061 ') for 
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Figure 6. The 68 and 95 per cent confidence contours in tlie wq ~ 
u>a plane determined from the /gas+CMB data (black, dashed 
contours) using the default dark energy model and 5 per cent 
systematic allowances. The solid red lines show the constraints 
obtained from the /gas data alone, using priors on Q^,h^, Q^^h^, 
Is, and h, as described in the text ("Section I5.5|l . The blue, dotted 
lines show the constraints from the /gas alone using priors on 
and h and assuming flatness. The figure shows how the 
CMB data contribute in constraining dark energy, especially at 
early times. 

other leading, planned (DETF Stage IV) dark energy exper- 
iments. We also find that, for the /gas experiment, the FoM 
can be boosted up by at least ~ 40 per cent by selecting an 
optimal redshift distribution of suitable clusters on which 
to carry out the /gas observations. Interestingly, the opti- 
mal redshift distribution of /gas measurments appears to be 
shifted towards low redshifts. 

As discussed in the text, a future /gas experiment will 
need to be preceded by a large X-ray or SZ cluster sur- 
vey that will find hot. X-ray luminous clusters out to high 
redshifts. A survey such as that planned with the Spectrum- 
RG/eROSITA mission should find several thousand of such 
clusters. Short 'snapshot' follow-up observations of the clus- 
ters with a new, large X-ray observatory should be able to 
identify a sample of ~ 500 suitable systems for /gas work. 
Attaining a precision of ~ 5 per cent with individual /gas 
measurements should be straightforward for an observatory 
with characteristics similar to Constellation-X, requiring ex- 
posure times of ~ 20ks on average. We note that the pop- 
ulation of galaxy clusters in the redshift, temperature and 
X-ray luminosity range of inter est has already been partially 
probed by the MACS survey (jEbeling et al.ll2001^ : Chan- 
dra observations of MACS clusters a r e used extensively in 
current / gas studies (|Allen et al.ll2004 iLaRoque et al.ll200^ : 
lAUen et al...2008. ) . The low- level of X-ray flux contamination 
from point sources observed in MACS clusters also alleviates 
the requirements on the instrumental PSF for dark energy 
work via the /gas method. 

In determining the predicted dark energy constraints. 



we have employed the same MCMC method used to an- 
alyze current data. The MCMC method encapsulates all 
of the relevant degeneracies between parameters and al- 
lows one to easily and efficiently incorporate priors and al- 
lowances in the analysis. We have included an array of such 
systematic allowances, with tolerances ranging from opti- 
mistic to pessimistic. Our technique differs from the DETF 
jAlbrecht et al.ll2006l ). who use a simpler Fisher matrix ap- 
proach in the prediction of dark energy constraints. Despite 
these differences, we have endeavored to make our calcula- 
tions of the FoM (Section 2) as comparable as possible. 

Benchmarking our results against those of the DETF for 
other, future 'Stage IV' dark energy experiments i.e. large, 
long-term missions, we find that the /gas experiment should 
provide a comparable FoM to future ground-based SNIa 
(FoM=8 - 22), space-based SNIa (FoM=19 - 27), ground- 
based BAG (FoM=5-55), space-based BAG (FoM=20-42) 
and space-based cluster counting (FoM=6— 39) experiments. 
Formally, the predicted FoM for the /gas experiment is 
comparable to 'pessim istic' scenarios f or we ak lensing ex- 
periments discussed bv lAlbrecht et al.1 (|2006l ). although the 
value falls short of the most optimistic DETF weak lensing 
predictions. The tight constraints on Qm and Qde for the 
/gas experiment will be of importance when used in combi- 
nation with other techniques. Interestingly, the 'pivot point' 
for the /gas experiment lies between those of the SNIa and 
BAG/weak lensing/cluster number count experiments, of- 
fering excellent redshift coverage in attempts to pin down 
the evolution of dark energy. 

We conclude that the /gas experiment offers a powerful 
approach for dark energy work, which should be competitive 
with and complementary to the best other planned dark 
energy experiments. 
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